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This paper deals exclusively with crossover designs for the pur- 
pose of comparing t test treatments with a control treatment when 
the number of periods is no larger than t + 1. Among other results 
it specifies sufficient conditions for a crossover design to be simul- 
taneously A-optimal and MV-optimal in a very large and appealing 
class of crossover designs. It is expected that these optimal designs 
are highly efficient in the entire class of crossover designs. Some com- 
putationally useful tools are given and used to build assorted small 
optimal and efficient crossover designs. The model robustness of these 
newly discovered crossover designs is discussed. 

1. Introduction. Crossover designs, where experimental units are used 
in two or more (p) periods for the purpose of evaluating and studying two 
or more (t) treatments, have proven effective in a wide range of applica- 
tions in agriculture, sensory testing and especially in phase I and phase II 
pharmaceutical clinical trials. The rigorous study of these designs and their 
optimality and efficiency has a history of at least 30 years. These optimality 
and efficiency results are almost all for those situations where all treatment 
comparisons are equally important. Some published works include Hedayat 
and Afsarinejad (1975, 1978), Cheng and Wu (1980), Kunert (1983, 1984), 
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Hedayat and Zhao (1990), Stufken (1991), Carriere and Reinsel (1993), Matthews 
(1994), Kushner (1997, 1998), Kunert and Stufken (2002) and Hedayat and Stufken 
(2003). Recently Hedayat and Yang (2003, 2004) have obtained additional 
results for balanced uniform designs and the designs suggested by Stufken 
(1991). We refer the readers to the excellent expository review paper by 
Stufken (1996) for additional references. 

In many pharmaceutical studies, experimenters are more interested in 
the comparisons between t experimental test treatments and an established 
standard or control treatment. In our study we shall designate the class 
of all such designs based on t test treatments, a control treatment and n 
experimental units each used in p periods by Qt+i,n,p- Unfortunately, the 
published literature contains very few useable results related to this very 
important problem. Perhaps the associated mathematical difficulties could 
be a major reason for this lack of useful published results. Majumdar (1988) 
considered A-optimal and MV-optimal crossover designs for comparing sev- 
eral test treatments with a control treatment and established that some 
known strongly balanced crossover designs can be combinatorially modified 
to obtain optimal designs for this problem. But Majumdar's result is limited 
to the situation where t = w 2 , p = m(w 2 + w) and m>2. The first nontrivial 
case will be t = 4 and p = 12. In crossover designs, for a variety of reasons, 
usually it is undesirable to have large numbers of periods. In many cases, 
the experimenter is interested in designs in which the number of periods is 
less than or equal to the number of total treatments, that is, p < t + 1. When 
p = 2, Hedayat and Zhao (1990) studied this problem and obtained useful 
results. 

Throughout this article, unless otherwise specified, we always assume 
3 <p < t + 1. The main purpose of this paper is to identify and construct 
crossover designs which are simultaneously A-optimal and MV-optimal in 
some desirable class of designs in terms of t test treatments versus a control 
treatment. Section 2 introduces the model and the notation. Section 3 con- 
tains four lemmas which provide useful mathematical tools for research on 
optimal crossover designs in general. The main result is presented in Sec- 
tion 4. Some assorted examples are given in Section 5. Section 6 discusses 
further issues and indicates some important open problems. The Appendix 
contains a series of propositions which are used in Sections 4 and 5. 

2. Model of response. Selection of an appropriate model of response for 
the data gathered under a crossover design in £lt+i,n,p is rather difficult and 
complicated. While several models have been introduced in the literature, 
the model which has been mostly entertained by design theorists is the 
traditional homoscedastic, additive and fixed effects model which, in the 
notation of Hedayat and Afsarinejad (1975), can be written as 

Ydku = A* + a k + 0u + T d(k,u) + Pd{k-l,u) + e ku, 
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(2.1) 

k = l,...,p; u = l,...,n, 

where Y dku denotes the response from unit u in period k to which treatment 
d(k,u) was assigned. In this model p is the general mean, at is the effect 
due to period k, (3 U is the effect due to unit u, T d rf. tU ) is the direct effect of 
d(k,u), Pd(k-i,u) is the- carryover or residual effect of treatment d(k — l,u) 
on the response observed on unit u in period k (by convention Pd(o,u) = 0) 
and the e^^'s are independently normally distributed errors with mean and 
variance a 2 . Recently another model has been introduced into the literature 
by Afsarinejad and Hedayat (2002) and some optimal designs have been 
identified. Kunert and Stufken (2002) used this new model of Afsarinejad 
and Hedayat and obtained additional optimal designs under this model. 

Hereafter we shall designate the t test treatments by l,2,...,t and the 
control treatment by 0. Throughout this paper, for each design d, we adopt 
the notation n diu , h diu , l dik , m dij , r di , r di and r d0 to denote the number 
of times that treatment i is assigned to unit u, the number of times this 
happens in the first p — 1 periods associated with u, the number of times 
treatment i is assigned to period k, the number of times treatment i is 
immediately preceded by treatment j, the total replication of treatment i in 
the n experimental units, the total replication of treatment i limited to the 
first p — 1 periods in the design and the total replication of control treatment 
limited to periods 2 to p, respectively. 

In matrix notation we can write model (2.1) for the np total observations 

as 

(2.2) Yd = fil + Pa + UP + T d T d + F dPd + e, 

where Y d = (Y dl i,Y d2 i,...,Y dpn y, a = (c*i, . . . , a p )', = (/3i, . . . ,(3 n )', r d = 
(r , . . . ,T t )', p d = (po,...,p t )', e = (eii,e 2 i,...,e pn ) / , P = ± n ®I p , U = I n ® 
t p , T d =(T^,...,T dn )' and F d = (F dv . . . , F' d J . Here T du stands for the 
p x (i + 1) period-treatment incidence matrix for unit u under design d and 
F du = LT du with the p x p matrix L defined as 

f 0ix(p-i) A 

\I(p-l)x(p-l) 0(p-l)xl/ 

The information matrix C d , for direct effects r d = (tq, . . . ,Tf)' , can now be 
expressed as 

C d = T dP r ± ([P\U\F d ])T d , 
where pr ± (Z) = I - pr(Jf) and pr(JT) = X(X'X)~X'. 
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3. Preliminary lemmas. For comparing test treatments with a control, 
the most frequently used optimality criterion is A-optimality. An A-optimal 
design minimizes 5Zj=i Var^fj — fo). Let Md be the information matrix for 
the contrast vector (t± — tq, . . . ,Tt — To)'. Then an A-optimal design for this 
contrast vector is a design which minimizes Tr(M7 ). Another optimality 
criterion, which is associated with A-optimality, is MV-optimality. An MV- 
optimal design minimizes Maxj=i j Var^(fj — fo). The following well-known 
lemma indicates that an A-optimal design is also an MV-optimal design. 

Lemma 1. An A-optimal design is also an MV-optimal design if its 
information matrix, Md, is a completely symmetric matrix. 

Lemma 2 points out a well-known algebraic relation between the two 
information matrices Cd and 

Lemma 2. The information matrix Md can be obtained from the infor- 
mation matrix Cd by 



Thus, Md can be simply obtained from Cd by deleting the first row and the 
first column of Cd- 

From Lemma 2, we can explore the algebraic properties of Md via Cd- 
But unfortunately, Cd for our problem is very complicated to deal with. It is 
known that when all test treatments and the control treatment are uniform 
in periods, that is, test treatments and the control treatment appear equally 
often in each period, we can ignore the impact of the period effects in Cd- 
This property is highly desirable when all test treatments and the control 
treatment are equally important. But when only the comparisons between 
test treatments and the control treatment are considered, this property is no 
longer desirable. In the latter case it is expected that the control treatment 
should have more replications than each test treatment. Lemma 3 shows 
that we can still ignore the period effects even if the replication of each 
test treatment and the control treatment are different. We state this result 
without a proof. 

Lemma 3. For any crossover design d, 



(3.2) C d < T' dW L {U)T d - T' dV T L {U)F d {F' dW L {U)F d )-F' dW L {U)Tdi 




with equality for a crossover design with Idik = r di/p, i 



= 0,...,t. 
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Notice that the uniformity in periods is just a special case of the stated 
condition. 

To find an A-optimal design, one needs to find a design that minimizes T^M^ -1 ). 
One way to achieve this is first to find the min^ g Q t+1 np Tf^Mj -1 ) and then 
characterize the design that achieves this minimum value. Although Lemma 2 
can help in simplifying the calculation of M^ 1 , it is still difficult to find the 
m i n den t +i „ P ^(M^ -1 ). The main difficulty is that for a general design d, 
one cannot write down the explicit expression of T^M^ -1 ). Lemma 4 shows 
a method for finding the explicit expression of the achievable lower bound 
for Tr (M7 1 ) for a very broad and general class of crossover designs. 



Lemma 4. For any design d in which < f^o < n (p — 1); we have 



where 



x 



t(np - r d0 - (l/p) g^i Eti n|J - (r d0 - (1/p) E"=i 

t(t-l) 



tp ( V ( m dii V n diu h di . 



r rf0 - f^doo + - >, n d0u n d o u 



xUt-1) 



n(p -l)(pt-t-l)-(pt-t + p- 2)r d0 +J2™ 2 d0u 



u=l 



yo = -[r do --Y^ 



n 



V 



dOu 



u=l 



^p[n(p - 1) - f d0 ] \m d00 n d0u h d0u j 

- (, P-1 if - " 

+ prdo r d0 r d0 - m d00 + - > n d0u n d0u 



x t 



np(p - l)f d o - r% - n(p h d0u 



Further, Tv{M^) 
satisfied: 



u=l 

+ ^ will hold when the following conditions are 
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(i) ldik = r d i/p, i = 0,...,t, 

(ii) T d pv (U)T d , T d pr ± (U)F d and F d pr- L (U)F d are invariant under any 
permutation of test treatments, 

(iii) each test treatment appears at most once in each of n first p — 1 
periods. 

Proof. It can be shown that the elements in T d pr ± (U)T c i at positions 
(i,i) and (i,j) (i / j) are r di - |ES=l n di« and -\YZ=l n diun d ju, respec- 
tively; the element in T' d pi L (U)F d at position (i,j) is m dij - ± Yl=\ ndiundju, 
the elements in F' d pi- L (U)F c [ at positions (i,i) and (i,j) (i ^ j) are r d i — 
\Yl=\ n \ iu and -\Yl=\ n diun d j u , respectively. 

Let S\ = It, S2, ■ ■ ■ , Sjsr be the set of all t X t permutation matrices on 
the set of t test treatments, where N = t\. Define M d = J2iLi S'iM d Si. By 
Lemma 2.2 of Majumdar and Notz (1983), M d is superior to M d under the 
A-criterion. Define 

Then Si = It+i,S2, ■ ■ ■ , Sn will constitute the set of all (t + 1) x {t + 1) per- 
mutation matrices on the set of all treatments leaving the control treatment 
unchanged. It can be easily verified that SlV'C d VSi = V Si C d SiV . 
By Lemma 2, we have 

1 N 1 N ~ ~ f 1 N 1 

(3.3) M d = -J2 S'iV'CdVSi = - J2 V'Si'dSiV = V' - E s/QS I V. 

i=l i=l I i=l J 

By (3.2) and Proposition 1 of Kunert and Martin (2000), we have 

TV TV 

E ^'O^ < E ^'(TdPr ± (^)7d 

i=l i=l 

-T l dV ^{U)F d {F' dV T ± (U)F d rF' dV ^{U)T d )Si 



(3.4) ^E^/^pr^^T,)^ 
j=i 

- (jz^{T' d V^{U)F d )s\ (Y^S'i(Kpr ± (U)F d )'S i 
\i=i I \i=i / 

x (^^\FdV^(U)T d )Sij. 

The equality sign in (3.4) will hold when (i) l di k = r di /p, i = 0, . . . ,t, and 
(ii) T^pr J -(J7)Trf, T d pr ± (U)F d and F d pi- L (U)F d are invariant under any per- 
mutation of test treatments. Then by utilizing the definition of Si, we observe 
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thatEiIi^'(^pr ± (^)r d )^, E J =i^ / (^pr ± (f/)F d )^ and (i^pr^tf) X 

Fd)Si have the form 

(3.5) ( a ™ f ™ Jlx * T ), m = 1,2,3, 

\Cm«Vjxl (o m — e m )ltxt + e m Jtxt / 
respectively. It can be shown that for YaLi Si {T' d pi L (U)T d )Si, 



ai = N\ r d0 

V P 



1 n \ 

^ w=l / 

&i = -r [np - r d0 - - ^ }^n diu , 
1 V ^«=ii=i / 



<* = /! = —, 

6i ai 
ei = -- — 7 + 



t-1 t(t — 1) : 



for EliSi'iT^^F^Si, 



a 2 = N \ m d00 V n d0u h dQu ] . 

. N*( 1» . 

t f-f \ p ^ 

1 = 1 \ ^ M=l 



02 

C2 = "T' 



t \ P V^ x 



e-2 



h + h 
t-l 



and for EiLi S- (F^^F^Si, 

a 3 =Nlr d0 - -J2n 2 d0u ), 



u=l / 

n t 



b3 = — [n(p-l)-f d0 --Y,Yj ) 

1 \ P u=li=l / 

Nf p-l_ 1 » _ 2 \ 
C3 = h = — ^do + - 2^ n rfo« . 
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n t 



iV M? - 1 _ 1 



*(*"!) I P P u =i 

~ ~ -{nip - 1) - f d o) + - 

P V i i 

r r U = l 1=1 



Thus, it is easy to see that when the control treatment is kept in the same 
positions, Ya=\ &i (^d 1 P T± (^)-Fd)Si will be maximized under the Loewner 
ordering when each test treatment appears at most once in each of n first 
p — 1 periods, that is, 

0.6) Y^liF'^mF<i%<{J Ch ~ e t hx \~ eJ ), 

r— f \ c 3^xi [03 — esjlfxt + esJtxt J 

w f t y, L — N[n{p-\)-f dQ ]{p-l) , - N (P-1 ? 1 y-n ~2 P~ 2 (V, v 

witn o 3 ^ ana e 3 - j(t±Tj\— r do p <L«=i n dOu ~\ nx 

(p— 1) — fdo)}. When < rao < n(p — 1), it can be easily verified that 03 > 0, 
63 > £3 and 0363 + (t — 1)0363 — £c§ > 0. So the right-hand side of (3.6) is a 
positive definite matrix and its inverse has the same matrix form as in (3.6) 
when 03, 63, c 3 and e 3 are replaced by 04, 64, C4 and e^. Here 



04 



63 + (t - l)e 3 



03^3 + (* - 1)0363 -*c| 



a 3 6 3 + (t - 2)a 3 e 3 - (t - l)c§ 



C4 



e 4 



^3 - e3)(a3&3 + (* - I)a3e3 - tcf) 
-C3 



0303 + (* - l)a 3 e 3 -tc\ 



c§ - 0363 



(63 - e 3 )(a 3 6 3 + (t - l)a 3 e 3 - ic§) 



The inverse matrix of the right-hand side of (3.6) can be expressed as D + 
C4J, where 



D 



(I4 — C4 0i x t 
fx i (h ~ e-i)hxt + (e4 - c 4 )J tX i 



Notice that C4 > due to c 3 < 0, and we have 

iV 



(3.7) x ^s/^pr^i^ ^/(^pr^T^ 



,i=l / \i=l 
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> ( Si{T' d w L {U)F d )i\ D ( £ sH(Fd prHU)T d )sh . 

The equality sign in the above inequality will hold when l dm = r d0 /p and 
each test treatment appears at most once in each of n first p — 1 periods. 
By (3.3), (3.4), (3.7) and by direct calculations, we have 

M d <M d < M d , 

and the equality signs hold when the three conditions in Lemma 4 are sat- 
isfied. Here M d = xl + y J, where 

x = —[bi-e 1 - (b 2 - e 2 ) 2 (6 4 - e 4 )] 

and 

y = —[e 1 - c 2 (a4 - c 4 ) - e 2 {b 2 - e 2 - / 2 )(&4 ~ e 4 ) - /|(e 4 - c 4 )]. 

Thus, to reach our conclusion, it is enough to show that the eigenvalues of 
M d are xq with multiplicity t — 1 and yo with multiplicity 1. 

Since M d = xl + yj, its eigenvalues are x with multiplicity t — 1 and 
x + ty with multiplicity 1. It is easy to check that 6 4 — e 4 = 1/(63 — £3), 
62 - e 2 = (t&2 + /2)/(* - 1) and 61 - ei = (t 2 6i - ai)/(t(t - 1)). By direct 
calculations, we obtain x = ceo- 

Now it suffices to show that x + ty = y$. Notice that b 2 + f 2 + (i — l)e 2 =0 
and by the definitions of a 4 , 6 4 , c 4 and e 4 , we have 

x + ty = —[h + (t - l)ei - tc 2 (o 4 - c 4 ) - t/|(e 4 - c 4 ) - /| (64 - e 4 )] 

_ oi tc|(6 3 + jt - l)e 3 + c 3 ) + /K03 + *c 3 ) 
iVi N[a 3 b 3 + (t- l)a 3 e 3 - tc 2 3 ] 

By direct calculations, we obtain 



0363 + (t- 1)0363 - tc. 



2 



TV* 
tp 2 



np{p - l)f d0 - rj - n{p — 1) 



f, 2 

n d0u 



u=l 



N 

b 3 + (t- l)e 3 + c 3 = —[nip - 1) - f d0 ] 

a 3 + tc 3 = . 

P 

Thus x + ty = y$. □ 
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4. Optimal crossover designs for test treatments versus a control treat- 
ment. In this section we will construct a family of crossover designs and 
prove their optimality for comparing t test treatments versus a control 
treatment over a large class of designs d in Q t +i,n,p for which Z^ofc = r do/p, 
k = 1, . . . ,p (i.e., the control treatment appears equally often in p periods) 
and nidii = for < i < t (i.e., no treatment is allowed to be preceded by it- 
self in d). The class of such designs is denoted by A. t +\,n,p- Before presenting 
the main theorem, we need some preliminary lemmas. To make our notation 
simple, we define the following four expressions: 




Ai = t(p-l)(np-r d o) -p[r d0 - 

(nt(p - 1) - tr d p - E"=i ndOundOu) 2 

n(p -l)( pt -t-l)-(pt-t + p- 2)fdD + Eu=i "k ' 



1^ 2 \ n (P~ 1 )(En= l n ^n d0 ^) 2 

ti nd ° U ) Mp-l)r d0 -rj - 
nt{p - 1) - tf d0 - E«=i ndoufidou 



Aa-P\ra np(p - l)r d0 - rj - n(p - l)E n u= A u 



n(p-l)(pt-t-l)-(pt-t+p-2)f d o + ES=i ™%u 
n i.P ~ !) Eu=i ndoufidou 



0, 



np(p - l)f d0 - rj - n(p - 1) £C=i n^ 0lt ' 



The next lemma provides a simpler lower bound than that in Lemma 4 
for any design d in A^ + x,n,p- This lower bound depends only on the property 
of the control treatment in the design d. 

Lemma 5. For any d £ At+i n p u>ii/i ./ixed r^o, we /iawe 
(4.1) ■&(M d - 1 )>* ( *- 1)2p + jg. 

xi = t(p- l)(np-r d0 ) -pfrdo - -min^n^ 0u J 

V P u=l ) 

[nt(p - 1) - tf d Q - minE"=i n d Q U ndQu] 2 

n{p -l){pt-t-l)-{pt-t + p- 2)f d0 + minY:™ =1 h 2 d0u 

and 

The equality in (4.1) ot/Z hold when the following two conditions are sat- 
isfied in addition to the three conditions in Lemma 4: 
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(i) n-diu, i = 1, ■ • ■ ,t, u = 1, . . . ,n, are either or 1, 

(ii) T,l=i n dOu> E«=i ndouhdou and E«=i«lo« are minimized. 

Here, the minimizations ofJ2u=i n d0w E"=i ndOufidOu andJ2u=i n d0u are 
over all possible designs d G A.fj r \^ n ^p with fixed r d0 • 

Proof. Since for any design d £ A t+ i np , with r^ = or np, the con- 
trasts between test treatments and the control treatment are not estimable, 
thus only designs with < r d o < n(p — 1) are considered here. When < 
r do < n Pi from Lemma 4, we know that Tr (jVf^~ ) > + with xq and 
yo as defined in Lemma 4. By Proposition A. 2 (for this and all other propo- 
sitions see the Appendix), we notice that tJ2i=iJ2u=i n diu n diu > t(n(p — 
1) - fdo) > J2u=i ndo u ndOu- We also know that m dii = 0, i = 0,...,t, and 
f d0 = ^^.Thn S 

1 

xq < -7- — tt-Ai, 
t(t-l)p 

1 . 

yo = — A 2 , 
tp 

and equality holds when n^, £ = 1, . . . , t, u = 1, . . . , n, are either or 1. So 
we have 

(4.2) Tr{Md -' ) > t ^l + f 2 

By Proposition A. 2, we have 

n 

n do u ndOu < nt(p - 1) - ifdo- 

For any d € At + i jniP with fixed r^o, let £i, £2 and £3 be independent variables 
satisfying the following inequalities: 



l d,0u 



£1 > nun ( 2 

min V n d0u n d0u < 6 < nt(p - 1) - tf d0 , 

By Proposition A.6, we know that min d (E"=i n d0u ), min d (E" = i n d0u h d0u ) 
and m.in d (J2u=i n d,0u) are fixed f° r given n, p and r^ - 
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Define Ai, A 2 , Oi and 02 to be the corresponding values of Ai, A2, 61 
and 6 2 , respectively, after replacing E2=i n dOw E«=i n d o u ndOu and E«=i ^dOu 
by £1, £2 and £3, respectively, in Ai, A 2 , 61 and 9 2 . Note that Eu=i n dOu> 
E™=i n d0und0u and E«=i ™dQu are n °t independent of each other, and thus 
the set {E"=i^do«' E«=i ™do«™do«, E«=i^ 0u } is just a subset of {£1, £ 2 , 
£3}. Thus for any d G Af +l niP with fixed r^ , we have 

f t(t-l) 2 P ,tp\ . (t^-Vfp t P 
mm t h -t— > mm h 

E: = x^.E: =1 "«o„,E: =1 n 2 do „v A ± A ^ A a a 2 

Define 

#(r<iO,£l, £2, 6) = ~ + ~~- 

Ai A 2 

To reach our conclusion, it is sufficient to show that for fixed r^o, H (rdo, £1, £ 2 , £3) 
is an increasing function of £1 , £ 2 and £3 . 

By direct calculation and Proposition A. 3, we have 

dtf(r d0 ,£i,6,6) _ t J} (t : i)V Q 



d& "VAi A? 

Also we have 

(43) M "-*^ M ^(|-^)> t 

5£ 2 VA| Af / 

The last inequality holds when ^ > due to Proposition A. 3 and (A. 10) 
in Proposition A. 5. Further, by applying Proposition A. 4 and (A. 9) in Propo- 
sition A. 5, the inequality still holds when ^ < jtj- 
By direct calculations, we have 

(44) 0H(r d0 M,^s) _ t J@ 2 2 (t-l) 2 ef 



96 VA2 A? 

And by the same argument, we have dH ( rd0 ^'&'& > q. Therefore for a given 



Ai 

minimized. Thus, the conclusion is obtained. □ 



r d0) ~ ^ ' £ + ^ will achieve the minimum value when £1, £ 2 and £3 are 
Ai A2 



Next we will introduce some definitions which are similar to the definitions 
in Kunert and Stufken (2002). A design d£ Q t +i,n,p is called: 

(i) a balanced test-control incomplete block design for the direct effects 
(with units as blocks) if (a) each test treatment 1 < i < t appears equally 
often in the design; (b) each test treatment appears in each unit at most 
once; (c) the number of units where any two test treatments % and j both 
appear is the same for every i 7^ j, 1 < i < t, 1 < j < t; (d) the control 
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treatment appears in each unit either {^] or [^] + 1 times; and (e) the 
number of units where the control treatment appears [^] times and test 
treatment i appears is the same for every 1 < i < t. 

(ii) a balanced test-control incomplete block design for the carryover ef- 
fects (with units as blocks) if the first p — 1 periods of the design form a 
balanced test-control block design for the direct effects in Qt+i,n.p-i- 

(iii) a balanced for test-control carryover effects design if (a) every test 
treatment 1 < i < t is immediately preceded by every other test treatment 
l<j<t equally often for every iy^j; (b) the control treatment is immedi- 
ately preceded by every test treatment 1 < i < t equally often; (c) every test 
treatment 1 < i < t is immediately preceded by the control treatment equally 
often, and no treatment including the control is immediately preceded by 
itself. 

(iv) a proportional frequency design for test-control on the periods if every 
test treatment 1 < i < t appears in every period exactly np ~j~p d0 times, and 
the control treatment appears in every period exactly ^ times. 

A design d* € Qt+i,n,p is called a totally balanced test-control incomplete 
crossover design if: 

(i) d* is a balanced test-control incomplete block design for the direct 
effects, 

(ii) d* is a balanced test-control incomplete block design for the carryover 
effects, 

(iii) d* is balanced for test-control carryover effects, 

(iv) d* is a proportional frequency design for test-control on the periods, 

(v) the number of units where test treatment j appears once in the first 
p — 1 periods and test treatment i appears once is the same for every pair 
i^ j; the number of units where the control treatment appears [— ] times 
in the first p — 1 periods and test treatment i appears once is the same for 
every test treatment 1 < i < t; the number of units where test treatment i 
appears once in the first p — 1 periods and the control treatment appears 
[^] times is the same for every test treatment 1 < i < t. 

Examples of such designs will be given in Section 5. The following lemma 
summarizes some useful properties of a totally balanced test-control incom- 
plete crossover design. 

Lemma 6. If d is a totally balanced test-control incomplete crossover 
design, then: 

(i) ldik = r d i/p, i = 0,...,t, 

(ii) Md is a completely symmetric matrix [since r^pr _L (?7)Trf ; T^pr _L ([/)i* 1 ( i 
and Fjpr- L (U)Fd are invariant under any permutation of test treatments], 
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(iii) En=i n do«> J2u=i n dOundOu and E"=l"do« are minimized over the 
designs with fixed r^ • 

We are now ready to summarize our findings in the following theorem. 

Theorem 1. For p < t + 1, a design d* is simultaneously A- optimal and 
MV-optimal in Aj + i injP if d* is a totally balanced test-control incomplete 
crossover design and rd*o minimize the right-hand side of (4.1). 

Proof. By Lemma 6, the conditions for the equality sign in Lemma 5 
hold. From Proposition A. 6, we notice that minY^" =1 n 2 d0u , min^" =1 n^n^ 
and minE"=i n do u are functions of rdo- Thus we have 



where x\ and y* are the corresponding x\ and y\ when d is d* . The last 
inequality holds due to Lemma 5. Thus design d* is A-optimal. And since 
M7, 1 is a completely symmetric matrix by Lemma 6, the MV-optimality of 
design d* follows by Lemma 1. □ 

5. Assorted mathematical tools useful for the construction of optimal 
crossover designs. While Theorem 1 specifies a set of sufficient conditions 
for a crossover design to be simultaneously A-optimal and MV-optimal, it 
will not provide any mathematical tools for constructing these designs. The 
purpose of this section is to alleviate this deficiency and present some math- 
ematical tools which in conjunction with some simple computer programs 
can help in constructing these designs. In practice, t and p are given and 
the job of the statistician is to find an optimal design for a given desirable 
n. If we want to rely on Theorem 1, then we are left with two tasks: first, 
to apply (A.13)-(A.15) in Proposition A. 6 into the right-hand side of (4.1) 
and run a simple computer program to find the rd*o which minimizes the 
right-hand side of (4.1), and second, to build a totally balanced test-control 
incomplete crossover design d where the replication of the control treatment 
in this design is r<i*o- 

Our experience shows that the family of totally balanced test-control in- 
complete crossover designs based on t, p and n contains very useful designs 
for our problems if r^o = n. We have observed that these designs are either 
optimal or at least highly efficient. Thus we shall concentrate our effort on 
this class of crossover designs. Not to leave the impression that the optimal 
designs are exclusively in this class of designs, we shall at the end of this 
section exhibit two optimal designs based on Theorem 1 with r^o ^ n. But, 
we should point out that the construction of optimal designs with r^ / n is 
not easy and we need more mathematical tools in this area. 
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Lemma 7. For given t, p, n and r^o = n, a totally balanced test-control 
incomplete crossover design exists only if both and are 

integers. 

Proof. Notice that the condition is an integer implies 2 is also 

an integer. This condition is necessary due to the uniformity for the test 
treatments on the periods. The condition that ^'~^^~^ n is an integer is 
also necessary since the design is balanced for test treatments. In fact, if 
the control treatment appears in the first or the last period, then any test 
treatment will be preceded by every other test treatment p — 2 times. And if 
the control treatment appears in any other periods, then any test treatment 
will be preceded by every other test treatment p — 3 times. So, totally we 
have [p— — 2)^ times that one test treatment is preceded by every other 

test treatment. Thus ^^rjy^ must be an integer. □ 

As we shall see soon, for many cases the necessary conditions in Lemma 7 
become sufficient for the existence of totally balanced test-control incomplete 
crossover designs. Lemmas 8 and 9 deal with cases for which p = t + 1 and 
p = t. As for the case of p < t, we shall provide a very useful tool that could 
help in developing such crossover designs. 

Lemma 8. For given t, p = t + 1, n and r^o = n, a totally balanced 
test-control incomplete crossover design exists if there is a balanced uniform 
design in &t+i,n,p- 

Proof. Simply replace treatment t + 1 in the balanced uniform design 
by the control treatment. Then it is easy to see that the modified design is 
a totally balanced test-control incomplete crossover design with the stated 
parameters. □ 

Notice that a necessary condition for the existence of a balanced uniform 
design in £lt+i n ,t+l is n = \(t + 1) for some positive integer A. According 
to Higham (1998), the class Qt+i,n,t+i contains a balanced uniform design 
when either n is an even multiple of t + 1 or t + 1 is a composite number, 
that is, it can be written as a product of two positive integers each larger 
than 1. 

Lemma 9. For given t, p = t, n and r^o = n, a totally balanced test- 
control incomplete crossover design exists if (i) S is an integer and t is a 
composite number, or (ii) is an even integer and t is a prime number. 
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Proof. We shall give an explicit way of constructing these designs. By 
the two conditions in the lemma, j is a multiple of t when t is a composite 
number and j is an even multiple of t when t is a prime number. Thus we can 
always construct a balanced uniform design in £lt,n/t,t- The totally balanced 
test-control incomplete crossover design with the stated parameters can be 
constructed in the following way: (i) Construct a balanced uniform design in 
£lt,n/t,t- (h) Replace treatment 1 by the control treatment in the balanced 
uniform design, (iii) Repeat step (ii) for the remaining test treatments to 
produce in total t new designs, (iv) Juxtapose the t new designs into a single 
design. The resulting design is the desired design. □ 

For the case p < t, our experience shows that the following steps lead 
successfully to a totally balanced test-control incomplete crossover design. 

Step 1. Construct a balanced incomplete block (BIB) design based on 
t test treatments in blocks of size p — 1. Let b be the number of blocks in 
this BIB design and number the blocks from 1 to b in an arbitrary fashion. 

Step 2. Construct p arrays each of size p by b. Fill the fcth array, 
k = 1, 2, . . . ,p, in the following way. Fill the entire b cells in the fcth row 
of this array by the control treatment. Note that each of the b columns of 
the array now has p — 1 empty cells. Fill these b columns arbitrarily with 
the b blocks of the BIB design. In this way we have produced p arrays each 
filled with the control treatment or the blocks of the BIB design. Juxtapose 
these p arrays into a big p by pb array. 

Step 3. Shuffle the positions of the test treatments in each column of 
the array produced in Step 2, for the purpose of converting the p by pb array 
into a totally balanced test-control incomplete crossover design. 

Although the preceding three steps cannot guarantee that in all cases a 
totally balanced test-control incomplete crossover design will be produced, 
our experience indicates that they are a very useful mathematical tool in 
building such designs when they exist. Indeed, we have succeeded in con- 
structing all small totally balanced test-control incomplete crossover designs 
or designs which are very close to such crossover designs. 

In the sequel we present some interesting examples based on the tools 
presented here. But first, a brief overview. For given t and p < t + 1 we 
can choose n so that the two integer conditions in Lemma 7 are satisfied. 
Then we can apply (A.13)-(A.15) in Proposition A. 6 into the right-hand side 
of (4.1) and run a simple computer program to find the r^*Q which minimizes 
the right-hand side of (4.1). If r^-o = n, we can utilize Lemmas 8 and 9 as 
the steps for the cases p < t and construct a totally balanced test-control 
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incomplete crossover design which is simultaneously A-optimal and MV- 
optimal by Theorem 1. In case r^* 7^ n i we can compare the corresponding 
minimum value of the right-hand side of (4.1) with the value when r^o = n. If 
the two values are very close, we can use Lemmas 8 and 9 as the steps for p < 
t to build a highly efficient crossover design. Ideally, we should investigate 
to see if it is possible to construct an optimal design when r^Q^n. 

We shall now concentrate on finding small size optimal designs for these 
practically desirable cases, namely (i) p = 3, (ii) p = 4 and (iii) p = 5. In 
addition, we exhibit two optimal crossover designs when r^*o ^ n. The gen- 
eral method of constructing optimal crossover designs when r^*o ^ n is very 
difficult and remains open. 

5.1. Simultaneous A-optimal and MV-optimal crossover designs for three 
periods. When t = 2, 3 or 4, then n must be a multiple of 3, 9 or 18, re- 
spectively. For the minimum values of n which satisfy the integer conditions 
in Lemma 7, we found that r^*o = n - Unfortunately, there are no totally 
balanced test-control incomplete crossover designs for either t = 2 and n = 3 
or t = 4 and n = 18. Thus, we searched for optimal designs for the next al- 
lowable value of n. We tried n = 6 for t = 2 and n = 36 for t = 4. Fortunately, 
r d*o = 6 and 36, respectively, for these situations and both designs can be 
easily constructed. For t = 3, the minimum n which satisfies the two integer 
conditions in Lemma 7 is 9. And indeed for this case r^*o = 9. Example 1 
exhibits one such optimal design for these parameters. 

Example 1. A-optimal and MV-optimal design for p = 3, t = 3 and 
n = 9: 

000231231 
123000123 
231123000 

When t = 5, n must be a multiple of 30 and indeed n = 30 satisfies the 
integer conditions in Lemma 7 and r^o = 30. Example 2 exhibits one such 
optimal design for these parameters. 

Example 2. A-optimal and MV-optimal design for p = 3, t = 5 and 
n = 30: 

000000000011453224351145322435 
231124535400000000002311245354 
114532243523112453540000000000 
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5.2. Simultaneous A-optimal and MV-optimal crossover designs for four 
periods. We tried to construct these designs for £ = 3, 4, 5, 6, 7, 8 and 9. 
The story is as follows. First, we note that the minimum n for these values 
of t must be a multiple of 4, 16, 40, 40, 28, 224 and 48, respectively. So we 
tried the minimum n. Fortunately, for all these cases r^*o = n. We did not 
attempt to construct the design for the case of t = 8 and n = 224 due to the 
size of n. For the remaining six cases we succeeded in constructing optimal 
crossover designs in the form of Theorem 1. For £ = 3 and n = 4, we can 
apply the technique in Lemma 8. When £ = 4 and n = 16, we can use the 
technique of Lemma 9. And for t = 5 and n = 40, t = 6 and n = 40, t = 7 
and n = 28 and t = 9 and n = 48, we successfully used the construction steps 
given after Lemma 9. We give samples of such designs in Examples 3-6. 

Example 3. A-optimal and MV-optimal design for p = 4, t = 5 and 
n = 40: 

0000000000135242531413524253141234543215 
1234543215000000000041352425315123454321 
5123454321413524253100000000002345132154 
2345132154352415314235241531420000000000 



Example 4. A-optimal and MV-optimal design for p = 4, i = 6 and 
n = 40: 

0000000000145462353135321612644266512345 
3532161264000000000042665123451454623531 
1454623531426651234500000000003532161264 
4266512345353216126414546235310000000000 



Example 5. A-optimal and MV-optimal design for p = 4, t = 7 and 
n = 28: 

0000000345671212345674567123 
12345670000000456 71233456712 
3456712456712300000001234567 
4567123123456734567120000000 



Example 6. A-optimal and MV-optimal design for p = 4, £ = 9 and 
n = 48: 

000000000000693174258456825963147123147528369789 
825963147123000000000000147528369789693174258456 
693174258456147528369789000000000000825963147123 
147528369789825963147123693174258456000000000000 
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5.3. Simultaneous A-optimal and MV-optimal crossover designs for five 
periods. We tried to build such designs for t = 4, 5, 6 and 7. For t = 4, n 
must be a multiple of 5 and the first n for which a totally balanced test- 
control incomplete crossover design with r^o = n exists is n = 10. For this 
value of n, we found that r^*o = n. The corresponding optimal design can 
be easily constructed using the tools in Lemma 8. However, we like to point 
out that if we need a design with bigger n for this case, then r^*o may not 
be n. In Example 10, we have exhibited an optimal design with n = 48 and 
r d * = 60. 

When t = 5, n must be a multiple of 25. So far we have not been able to 
construct a totally balanced test-control incomplete crossover design with 
r do = n when n = 25. However, we can construct such a design when n = 
50 by using the tools in Lemma 8. Notice that when n = 50, rd*o = 60 
minimizes the right-hand side of (4.1) and the corresponding minimum value 
is 0.24179. Although r^o = 50 does not minimize the right-hand side of (4.1), 
its corresponding value is 0.24299, which is 99.5% efficient relative to the 
minimum value. So this design is highly efficient or even optimal. 

When t = 6, a totally balanced test-control incomplete crossover design 
will be relatively very large since n must be multiple of 75. We tried n = 
30 knowing that we cannot use Theorem 1 to conclude optimality, but we 
hoped for a very efficient design. For n = 30, we found rd*o = 30 and the 
corresponding minimum value of the right-hand side of (4.1) is 0.55044. 
We used the construction steps after Lemma 9 and found a highly efficient 
design. This design is given in Example 7. For this design d, its Tr(M7 ) = 
0.55419, which is 99.3% efficient relative to the minimum value. So this 
design is highly efficient or even optimal. 

Example 7. Efficient design for p = 5, t = 6 and n = 30 with efficien- 
cy =99.3%: 

000000612345234561345612561234 
123456000000612345561234612345 
612345234561000000234561234561 
234561561234345612000000123456 
561234345612456123612345000000 

When t = 7,n must be a multiple of 35. So far we have not succeeded in 
constructing an optimal design in the form of a totally balanced test-control 
incomplete crossover design with r^o = n = 35. However, we have been able 
to construct such a design for n = 70. For n = 70, we found r^» to be 70 and 
we shall present one such design in Example 8. The design in Example 8 is 
split into two parts d\ and d<i each based on n = 35. While the union of the 
two designs is optimal, each of them is a highly efficient design in the class 
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of designs with n = 35. When n = 35, the minimum value of the right-hand 
side of (4.1) is 0.61529 while T^M^ 1 ) = 0.61904 and Ti(Mj^) = 0.61927. 
Thus d\ is 99.39% efficient and cfo is 99.36% efficient. 

Example 8. A-optimal and MV-optimal design for p = 5, t = 7 and 
n = 70: 

d = d\ U d% , 

where 


2 3 4 5 
di = 6 7 1 2 
5 6 7 1 
7 12 3 

and 
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Before closing this section we shall present two totally balanced test- 
control incomplete crossover designs for which r^*o 7^ n - Such designs are 
extremely difficult to construct. For p = 3, t = 7 and n = 49, we have rd*o = 
42 and for p = 5, i = 4 and n = 48, we have r^*o = 60. 

Example 9. A-optimal and MV-optimal design for p = 3, i = 7 and 
n = 49: 

0000000000000012345674567123456712312345674567123 
3456712123456700000000000000123456734567123456712 
4567123456712334567123456712000000000000001234567 

Example 10. A-optimal and MV-optimal design for p = 5, f = 4 and 
n = 48: 

000012341234123400001234123412341234000012341234 
123400002413241312340000241324132413241300002413 
241324130000314224132413000031423142314231420000 
314231423142000031423142314200000000432143214321 
432143214321432143214321432143214321000000000000 
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6. Discussion and closing remarks. In this article we imposed two con- 
ditions on the class of competing designs: (i) In each design the control 
treatment should appear equally often in all p periods, (ii) In each design 
no treatment is allowed to be immediately followed by itself in any exper- 
imental unit (i.e., ma = 0, i = 0, 1, . . . ,t). Although the optimality of the 
newly discovered designs is over this restricted class, we are highly confident 
that these designs are very efficient over the entire class of designs. Indeed, 
for some t, p and n these designs could be optimal over the entire class. 
To support this optimism, we can see from Lemma 3 that for any design d 
in which condition (i) is not satisfied, (3.2) may become a strict inequality. 
It is possible that this condition is a necessary condition for the optimal 
design in the entire class and therefore the restricted condition (i) could be 
removed. As for condition (ii), if nidu > for some 1 < i < t, xq in Lemma 4 
may not yield a significantly bigger value since the gain in the second part 
will be reduced by the loss in the first part. Consequently, we think the 
restricted condition man = for all 1 < i < t could indeed be a necessary 
condition for the optimal design, or at least the gain with the condition 
removed is very little. Further, if m^oo > 0, the numerator of 02 becomes 
n(p — l)(J2 n d0uf l d0u — P m doo) an d @2 could be zero or negative. Thus, the 
corresponding inequalities (4.3) and (4.4) may not be greater than zero for 
a general design d and consequently Ti^M^ 1 ) may not be minimized when 
J2 nd0v,nd0u or ^dOu i s minimized for fixed rdo ■ In this situation, we may 
have to consider the relationship among J2 n do u > S n rfo«"-dOit and ^dOu ■ We 
may find optimal designs for some special parameters, but we feel it will be 
very difficult to find optimal designs for the general parameters p < t + 1 . 

It is natural to postulate that there could be better designs than those we 
have identified in this paper. While we do not have any general evidence for 
that, what so far we can say is this. The designs characterized in this paper 
are highly efficient if not optimal. We have substantial numerical evidence 
in support of this. Here is a typical example from the assorted examples 
that we have produced in our ongoing research in this area. Consider the 
case t = 7, p = 4 and n = 28. We used a computer along with some algebraic 
methods and searched for the possible lower bound of Tv{Mj l ) under the 
subclass of designs in which the control treatment appears equally often 
in p periods. We found that there could be a design with Ti{MJ l ) equal 
to 1.02252. While we are not sure if there is a design with such a trace, 
let us assume that there is one. Clearly this hypothetical design is better 
than the A-optimal design d* with Tv{Mj l ) = 1.02327 which we displayed 
in Example 5 within the subclass ht+i )n)P of designs for these t, p and n. 
But note that d* in Example 5 is at least 99.9% efficient in the larger class 
without the restriction ma = 0, i = 0, 1, . . . ,t. This is not an isolated case, 
and as we mentioned, we have observed this phenomenon in many cases. 
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Another important issue worth discussing here is the status of the model 
robustness of the optimal and efficient designs discovered in this article. We 
used carryover Model (2.1) while searching for optimal and efficient designs. 
It is quite possible that upon data collection and analysis we may discover 
that another model could be more appropriate than Model (2.1). Conse- 
quently, the optimal design which was used under the postulated Model (2.1) 
might no longer be optimal or even efficient under the model specified in light 
of the data. Thus it is prudent to recommend a design to the experimenter 
which is optimal or at least efficient under several different likely models. For 
crossover studies several simple and lower models than Model (2.1) have been 
used in practice. We shall explore here how our designs in this article perform 
under the model without carryover effects (two-way elimination model), the 
model which contains only direct treatment effects and subject effects (one- 
way elimination model) and the model which contains only direct treatment 
effects (zero-elimination model). In Table 1 we have listed 15 designs of 
which 10 have already been displayed in Section 5 and the remaining 5 can 
be obtained by the procedures in this paper. This table lists the efficiency 
of each design under the carryover Model (2.1) in the subclass Af+i ;nj p of 
designs (e c ), as well as zero-way (eo), one-way (ei), and two-way (e2) elim- 
ination models in the unrestricted class of designs. The efficiencies eo, e± 
and e2 are based on the result of Hedayat, Jacroux and Majumdar (1988). 
It is clear that for most of these designs, if the design is A-optimal under the 
carryover model, then it is also A-optimal under one-way and two-way elim- 
ination models, and highly efficient under the zero-way elimination model. 



Table 1 

Efficiencies under different models 



Design 


P 


t 


n 


r 


e c (%) 


eo (%) 


ei (%) 


e 2 (%) 
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100 


97.50 


100 


100 
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100 


100 


100 


100 
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30 


36 


100 


100 


100 


100 
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30 


30 


100 


99.84 


99.85 


99.85 
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49 


42 


100 


99.98 


99.97 


99.97 
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100 


94.4 


98.75 


98.75 
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16 


16 


100 


96.55 


100 


100 


8 
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40 


40 


100 


98.18 


100 


100 
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4 
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40 


40 


100 


99.16 


100 


100 
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28 


28 


100 


99.82 


100 


100 
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48 


48 


100 


100 


100 


100 


12 
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48 


60 


100 


96.43 


97.56 


97.56 


13 
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50 


50 


99.50 


93.10 


96.87 


96.87 


14 
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30 


30 


99.30 


95.23 


98.53 


98.53 


15 
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7 


70 


70 


100 


96.68 


99.47 


99.47 
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Even an efficient design under the carryover model remains highly efficient 
under the other three models. Notice that while Design 13 is highly efficient 
under the carryover model, it is not that efficient under other models. The 
reason is this: For computational simplicity this design was constructed un- 
der r^o = 50 while our computation showed that r^*o is 60. Therefore, it 
is very safe to conclude that optimal and efficient crossover designs which 
are constructed in this article remain optimal or highly efficient under lower 
case models which are discussed here. Our research effort in this area is 
continuing. 

APPENDIX 

PROPOSITION A.l. For any deA t+ltriiP , we have 

nip — 1) 

Proof. Since the control treatment appears equally often in periods, 
we have f^o = (p — l)^oi- Also Z^oi < n/2 since no treatment (either test 
treatment or control) is followed by itself. Thus we obtain the conclusion. 
□ 



Proposition A. 2. For any d € A t+ i jniP , where p<t+l, we have 

n 

n do u ndOu < t[n(p - 1) - fdo]- 



u=l 



Proof. For given r d o, since £™=i n d o u = fdo and < n d0 « < |, we have 

(A.i) E4>«< X' whenp - 4 ' 

u=i lf,jo, whenp = 3. 

On the other hand, we have J2u=i n d0v,nd0u < E«=i ™do«(™dOu + !)• Notice 
that 3 < p < i + 1; then by applying (A.l) and Proposition A.l, we obtain 
the conclusion. □ 

In the next three propositions, £i, £2, £3, Ai, A2, ©1 and ©2 have the same 
definitions as those in the proof of Lemma 5. Notice that for any d £ A.t+\, n ,p 
and given r^o, we have 



6 >nun( J2 n d0u) > r do, 

\ u=l / 

nt(p-l)- if do > 6 > mjn ^ n do u ndo u >?do, 
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£ 3 >irun( Y.™dOu) 



Proposition A. 3. For any d € A t +i, njP , where 3 <p < t + 1, we have 
(A.2) Ai>(t-1)A 2 . 

Proof. Since nt(p — 1) — tf d Q > £2 > ^do and £3 > f d o, so we have 

(nt{p-l)-tr d0 -£ 2 ) 2 

n(p - l)(pt - 1 - 1) - {pi - 1 + p - 2)r d0 + £3 

(nt(p - 1) - (t + l)rdo) 2 



< 



and 



n(p - l)(p* -t-l)-(pt-t+p- 3)f d0 



n(p - l)(t - 1)(6) 2 > n{p - l)(t - l)fdo 



np(p - l)r d0 - - n(p - 1)6 "(p - I) 2 - r d0 

2 

We also notice that tp(r d0 — < t(pr d Q — -Jf)- By direct calculation it is 
sufficient to show 

r do\ -,\, \ n(p- 1)^-1)^0 



* (p^do - ^) - *(P - l)(np - r d0 ) 



n(p - l) 2 - r d0 

(A.3) 

, (nt(p-l)-(t + l)r d0 ) 2 <Q 

n(p - l)(pt - t - 1) - (pt - t + p - 3)r d0 ~ 

The left-hand side of the preceding expression can be written as 

-n Up - ^ )( P - 1 - r J») + {p - 1){t - l){ ^ ,n) 



n J V n J (p- l) 2 - (fdo/n) 

(t(p-l)-(t + l)(f d0 /n)) 2 

(p - l)(pt - t - 1) - (pt - t +p - 3){f d0 /n) J • 

From the proof of Proposition A.l, we know that ^ < ^ + 1/2. Define 
x = — ; then ^ < x + 1/2. Simple algebra can show that (A.3) is equivalent 
to f(x) > when < x < (p — l)/2, where 

/(x) = i(2p - 1 - 2x)(2p - 3 - 2x) + 4(t ~ 1)x 

P- 1 

4(t(p-l)-(t + l)x) 2 



(p- l)(pt-t- 1) - (pt -t + p-3)x 
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It can be shown that (p — 1) {pt — t — 1) — (pt — t + p — 3)x > when < x < 
{p — l)/2, so it is equivalent to showing g(x) > when < x < (p — l)/2. 
Here 



9{x) 



t(2p - 1 - 2x)(2p - 3 - 2x) + 



4(t- l)a 



x [(p- -t-l)- (pt-t + p-3)x\ -A(t(p- 1) - (t + l)x) . 



It can be checked that g(0) > and g{^-) > 0. If we can show that g(x) is 
a monotone function when < x < (p — l)/2, then we have reached the 
conclusion. By direct calculation, we have 



g'(x) = -12t(pt -t+p- 3)x 2 



+ 



t 2 (2Ap 2 - A8p + 8) + t(16p 2 - 72p + 40) + 



16(* — 1) 



p-1 



- t 2 (12p 3 - 36p 2 + 23p + 1) - t(Ap 3 - 28p 2 + 39p - 9) + 4. 
Notice that g'{x) is an increasing function when x G (— oo,M), where 
t 2 {2Ap 2 - 48p + 8) + t(\Qp 2 - 72p + 40) + 16(£ - l)/(p - 1) 



M 



2At(pt - t + p - 3) 



It can be verified that M > (p — l)/2. So <?'(:e) is an increasing function when 
x G [0, (p — l)/2]. Also we can verify that 



t 2 (-3p 3 + 9p 2 -Ap-2) + t(p 3 - p 2 - Ap + 6) - 4. 



Notice that -3p 3 + 9p 2 - Ap - 2 < when p > 3; thus 

c/' (^y^) < 2i(-3p 3 + 9p 2 -Ap-2) + t(p 3 -p 2 -Ap + 6)-A 

= t(-5p 3 + 17p 2 - 12p + 2) - 4 
<0. 

So we have </(x) < when x G [0, (p — 1)/2]; thus ^(x) is a monotone function 
when x6 [0,(p-l)/2]. □ 

Proposition A. 4. For anyde At+i,n, P , withp< t+1, and ^ G [0, f=£], 
we have 



(A.4) 



Ai 



> 



*(p-l) 



(t-i)A 2 - t(p-i)-r 
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Proof. Since nt{p — 1) — tr d0 > £2, Ai will be minimized when £1, £2 
and £3 are minimized, that is, 



(A.5) 



Ai > t{p - l)(np - r d0 ) -ip- l)r d0 

[nt(p-l)-(t + l)f d0 ] 2 



nip - l)(pt-t- 1) - ipt -t + p- 3)fdo 



Also we can see that A2 will be maximized when £1, £2 and £3 are minimized, 
that is, 



(A.6) 



A 2 <(p-l)r d0 



nip - l)r d0 
nip - l) 2 - fdo 



It suffices to show that 



1 



n 



[{t(p - 1) - l)Ai - tit - l)(p - 1)A 2 ] > 0. 



Applying (A.5) and (A.6), and noticing that r d o = r^jf^o, we can show that 



-[(t(p-l)-l)Ai-t(t-l)(p-l)A 2 ] 



>[t(p-l)-l] 



(A.7) 



pt(p-l)-p(t + l) 



r d o 



n 



t 2 {p~l) 



{p-l){pt-t-l)-{pt-t + p- 3){r d0 /n) 



tp{t-l){p-l) r -^. 

n 



, and the right-hand side of (A.7) is a decreasing function 



Since < f -f < B=± 
of fdo, we further have 

1 



(Ai 



[(*(p-l)-l)Ai-t(t-l)(p-l)A 2 ] 
n 

>[t(p-l)-l] 

pi(p - 1) - p(p - 1) 

i 2 (p-l) 



(p-l)(pi-t-l)-(pt-i + p-3)(p-l)/(t + l) 
p — 1 



ip(i-l)(p-l) 



t + 1 



When p > 4, by applying the condition that p < t + 1, we can show that 
(pt - 1 - 1) - pt ~*+f~ 3 > , t (p - 1) - 1 > and p(p - 1) {t - 1) - 3t > 
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p(p ~ 1 3 )(t ~ 1) , so by (A.8) we have 

-[(t(p-l)-l)Ai-i(t-l)(p-l)A 2 ] 
n 



> [t(p - 1) - l][p(t - - 1) - 3f] - tp(t - l)(p - 1)5-1 



2*o.-i)'(«-»(i-rri) 



>0. 

When p = 3 and t > 3, notice that 1 + ^ < f. By (A.8) we have 

-[(i(p-l)-l)Ai-t(t-l)(p-l)A 2 ] 
n 

x / , 2i 2 t(t-l) 
> 6(2* — l)(t — 1) — (2t - 1) 12- v 



2t - f t + l 
> (2t-l)(4t-6)-3t(t-l) 
>0. 

When p = 3 and t = 2, due to (A. 5) and (A. 6), we have Ai > 12n — 6r d o 
(4 f T~ 3 f-° )2 and A 2 < 2r d0 - ^4^. Thus we have 

-[(t(p-l)-l)Ai-t(t-l)(p-l)A 2 ] 
n 

> 36 _ 26 r^ _ 3 (4 ~ 3 ^/^) 2 + 8r */" 



n 6-4f d0 /n i-r d0 /n 

>36-37^-3 (4 " 3rWn)2 
n 6-4r d0 /n 

>0. 

The last inequality can be easily verified when < ^ < | . □ 
PROPOSITION A. 5. For any d e A t+ i |n)P; mi/i p < t + 1, we have 

(A.9) o<^ §1 <e, 

Furthermore, when — > ttt, we have 

(A.10) e 2 >9i. 
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Proof. It is easy to verify that 0i > 0. We only focus on the remaining 
inequalities. By an argument similar to that in Proposition A. 3, we have 

nt(p - 1) - tr d0 - & 



n{p - l)(pt -t-l)-(pt-t + p- 2)r d0 + £ 3 



< 



nt(p - 1) - (t + l)f d0 



n(p - l)(pt - t - 1) - (pt - 1 + p - 3)f rf0 



and 



n(p - 1)6 



> 



n{p — 1) 



np(p - l)f d0 - rj - n(p - 1)£ 3 n(p - l) 2 - f d0 ' 
As for (A. 9), it is sufficient to show that 

p-1 (p- l)(pt -t-1)- (pt-t + p-3)r d0 /n 



(A.ll) 



(p - l) 2 - f d0 /n 



t(p-l)-(t + l)f d0 /n 



> 



t(p 



t(p-l) 

Direct calculations show that (A.ll) is equivalent to 



(t + l)[t(p -1)-1] 



,r d0 



[( p -iy^+t-i)-t(p-i)Y-^<o, 

n 



which holds when < % < (p ~ffi ~m 2 +t~i)-t] 
„ ( P -i){(p~i)(t 2 t] . d 



( t+i)|t(p-i)-ij Notice that {p-l){t 2 +t- 



l)-t>t 2 (p-l)-t, so ""(^gl'grj^^ > From Proposition A.l, 

we know that %^ < 2=1; thus (A.9) holds. 
For (A. 10) it suffices to show that 

t(p-l)-(t + l)f d0 /n 



(A.12) 



p-1 



> 



(p- l) 2 -f d0 /n (p - l)(pt - t-1) - (pt- t+p-3)r d0 /n' 
Direct calculations show that (A.12) is equivalent to 



1)(^\ 
n J 



(t + 2)(p-l) r -^ + (p-l) 2 <0, 
n 



which holds when ^ £ [2— j- iP ~ F rom Proposition A.l we know that 
< 2=1. Thus (A. 10) holds when ^ > 2=1. □ 

Proposition A. 6. For any design d G A( + i injP and given r d Q, 
(A. 13) mjn ^ ^ n 2 d0 ^j = r d0 + (2r d0 - n) 



>do" 




'r d Q 




— n 




n 




n 
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(A.14) 



and 
(A.15) 



min ^2n d0u n d0u 



,u=l 



rao + 2f d0 - n + 



p-1 



n 



2f d0 + 



p-1 



when 

n + ( 2r d0 -2n + 



- n 
fdO 

p — 1 

fdO 



rdo 



n 



<n-f d0 + n 



rdo 



n 



p — 1 

otherwise, 



n 



n 



rdo 



n 



mm 

d 



V 11=1 



"^ ~ 




~fdo~ 




— n 




n 




n 



J2u=i n dOu> E«=i ndoufi-dou and J2u=i n dOu can achieve their minimum values 
when d is a totally balanced test-control incomplete crossover design. 



Proof. It is straightforward to show (A. 13) and (A.15). Here we will 
prove (A.14) only. First we notice that 

n n 
^ ndOufldOu = J! ™dOu + E] ™dOu- 

u=i u=i «gr 

Here T is the set of units which receives the control treatment in the last 
period. So there are r d o — fdo = 3; units in T. 

For any hdo u , u = l,...,n, 2~2u=i n dOv,n d Q U will be minimized when we put 
the ^jjfj- smallest values among all n d ou into T. Next, we will show that for a 
given fdo, 2~2u=i n dQvJ l dau will be minimized when d is a balanced test-control 
incomplete block design for the direct effects and carryover effects. 

Suppose that there are some flow's, say hdoi and hd02, such that n d oi — 
nd02 > 2. Then we can replace h d oi by n' d01 = n dm - 1 and h d Q 2 by h' d02 = 
n d Q2 + 1 and keep the others unchanged. Direct calculations show that the 
value of 2^u=i™dOu * s decreased by at least 2. Meanwhile, Ylu&v^dOu is in- 
creased by at most 1. So J2u=i n d0und0u is decreased by at least 1. Thus 
2~2u=i n d0und0u w ih be minimized when d is a balanced test-control incom- 
plete block design for carryover effects and T is the set of units which has 
the smallest values among fidOu- When d is a balanced test-control in- 
complete block design for the direct effects and carryover effects, it satisfies 
this condition. By direct calculations, we obtain the conclusion. □ 
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